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Abstract 

We give the conversion equations which lead from experimental values of the 3p 
Is transition energy in pionic hydrogen and the total width of the Is level to values 
of the s-wave threshold scattering parameters for the processes iT~p n^p and 
7r~p TT^n respectively. Using a three-channel potential model, we then calculate 
the electromagnetic corrections to these quantities, which remove the effects of 
the Coulomb interaction, the external mass differences and the presence of the 
7n channel. We give the s-wave scattering parameters obtained from the present 
experimental data and these electromagnetic corrections. Finally we discuss the 
implications for isospin invariance. 
PACS: 13.40.Ks; 13.75.Gx; 36.10.Gv 
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1 Introduction 



Accurate values of the s-wave scattering parameters Qcc and Oco, for tlie pro- 
cesses Tx^p —>■ Ti^p and TT^p —>■ Ti^n at tlie n^p tiireshold, can be obtained 
from measurements of an np Is transition energy and the width of the Is 
level in pionic hydrogen (H) via a refinement of the DGBT formulae [I]. As 
the most accurate experimental data available are for the 3p — *■ Is transition 
energy, we will consider only this case in the present paper. The most recent 
fully published experimental values of these quantities are given in Ref. |2|. A 
successor experiment [3] is nearing completion, and preliminary results from 
this experiment for the transition energy [3] and for the width [5] are available. 
In Section 2 of this paper we shall give a detailed derivation of the relevant 
conversion equations which connect the experimental data with the scatter- 
ing parameters, paying particular attention to the sources of uncertainty and 
making clear which pion-nucleon {ttN) interactions contribute to the param- 
eters. We use these conversion equations to obtain the values of Occ and Oco 
from the results in Refs. [^115] . These equations can also be used for the final 
results of the experiment [5] and for those of any future experiment. 

The conversion equations just mentioned are based on the formalism of Ref. [6]. 
It was the first work to generalise the systematic derivation of the DGBT for- 
mulae for a one-channel system [7118] to a multichannel system, and to apply 
the results to pionic H. This systematic derivation is based on the general 
principles of scattering theory and on the analytic continuation of the matrix 
K (defined later in Eq. ([3])) into the complex s-plane, to those (complex) val- 
ues of s which correspond to the levels of pionic H. All this is independent of 
any phenomenological model of the dynamics of the system. The formalism 
of Ref. [0] enables one to extract directly and unambiguously from the exper- 
imental data the values of precisely defined threshold scattering parameters 
and ttco, whose errors are dominated by the experimental errors of these 
atomic data. Then, as a second quite distinct step, a dynamical calculation of 
the electromagnetic (EM) corrections to these quantities is needed if we are 
to obtain hadronic scattering parameters. 

In Section 3 we shall describe our calculation of the electromagnetic (EM) 
corrections to the threshold parameters, which uses a three-channel potential 
model to remove the effects of the Coulomb interaction, the external mass 
differences and the presence of the jn channel. We shall also give the values, 
together with estimated errors, of the corrected scattering parameters (de- 
noted by Sec and a^o) which result from and a^o- In Section 4 we shall 
set our calculation of the EM corrections in the context of previously pub- 
lished calculations, and attempt to summarise our present understanding of 
these corrections. We shall see that, while the experimental data from pionic 
H yield extremely accurate values of ttcc and ttco, the calculation of the EM 
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corrections needed to obtain hadronic quantities is subject to a great deal of 
uncertainty. 

In Section 5 we shall use the value of the s-wave scattering length a+ for 
TT+p elastic scattering, obtained from a recent phase-shift analysis (PSA) [H] 
of the data up to a laboratory pion kinetic energy (T^) of 100 MeV, which 
employed the EM corrections of Ref. [10], also calculated using a potential 
model. To ensure that the two calculations of EM corrections were consistent, 
we adopted, for the hadronic potentials in the calculation described in Section 
3, the same forms as were used for the corrections given in Ref. [TU]. We shall 
show that this value of a+ shows a large and significant difference from the 
value of dec + V^dco, calculated from the most recent values of dec and dco 
given in Section 3. 

This reinforces the evidence for the violation of isospin invariance given in 
Ref. [9], from the inability of the phase shifts, obtained from the analysis of ir^p 
elastic scattering data alone, to reproduce the data on n'p charge exchange 
scattering. As discussed in Ref. [H], at least part of this violation of isospin 
invariance has to be attributed to the presence of residual EM effects in the 
phase shifts obtained from the experimental data after the application of our 
corrections. This is because the EM corrections calculated using a potential 
model are only partial (stage 1) corrections; further (stage 2) corrections are 
needed to remove these residual EM effects. 

We do not include in the present paper the information from the accurate 
measurement of a transition energy in pionic deuterium (D) [TT]. There is an 
extensive literature on this subject [T2|13fl4j . The most recent of these papers 
[H] concludes tentatively that, when stage 2 corrections (calculated using chi- 
ral perturbation theory at leading order) are taken into account, the atomic 
data from pionic H and pionic D taken by themselves satisfy isospin invariance 
at the hadronic level. However, Ref. [11] emphasises that an improved mea- 
surement of a transition energy in pionic D is needed, as is a major effort on 
the theoretical side, to extract a scattering parameter with better precision. 
This complex subject lies outside the scope of this paper. Our aim here is to 
present what can be learned directly from the experimental data for pionic H 
and to relate the results consistently to those obtained from our PSA [9J. 



2 Conversion formulae 

The first piece of experimental data from pionic H is the 3p Is transition 
energy W{3p Is). The positions of the 3p and Is levels relative to the ir'p 
threshold we call W{3p) and W^(ls) respectively. For the numerical calcula- 
tions we will use always the most recent values of the physical constants [T5]; 
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errors will be given only if they are relevant for our purpose. 

Nonrelativistically the point charge Coulomb potential V^'^ contributes 
—359.438 eV to W{3p). The shift due to vacuum polarisation is —0.011 eV 
[TB] . Furthermore the 3p level is split by the hyperfine interaction and the 
relativistic effects. The former is discussed on page 413 of Ref. [6], while the 
latter is given by Eqs. (9,10) of Ref. [17]. For j = 3/2 these two effects add 
to —0.0013 eV, for j = 1/2 to —0.0055 eV. With the present experimental 
resolution of W{3p —>■ Is) this splitting cannot be seen. We therefore take 
a weighted average of —0.003 eV, with an uncertainty of 0.002 eV, as the 
contribution of these effects to W{3p). The shift in the 3p level due to the 
hadronic interaction, including the effect of the extended charge distributions 
of vr^ and p, is negligible. Thus 

W{3p) = -359.452(2) eV. 

Since in the transition 3p Is the kinetic recoil energy of the tt" atom in the 
cm. frame is 0.004 eV, we have 

W{ls) = -W{3p Is) - 359.456(2) eV . (1) 

Nonrelativistically V^"^ contributes —3234.945 eV to W{ls). The relativistic 
shift for an infinitely heavy proton is —0.215 eV. This comes either from 
Ref. [T7j or from integration of the Klein-Gordon equation [16j. In addition, 
both these references give a relativistic recoil shift of +0.037 eV and a con- 
tribution of +0.010 eV due to the magnetic moment of the proton. Ref. [16] 
also quotes a contribution of —0.018 eV for vacuum polarisation at 0{a^) 
and a vertex correction of +0.007 eV. The potential V^^ for vacuum polarisa- 
tion at 0{a^) is the Uehling potential, modified at short distances due to the 
finite extension of the charge distributions of tt" and p. Adding V"^ to V"^^, 
our calculations give a shift of —3.244 eV. The finite extension of the charge 
distributions also modifies V'^'^ itself through the addition of a potential 

2 

y^^* = -a/r(er/(r/c) - 1), with = -(< >^ + < >p) . 

3 

The rms charge radii of vr" and p were taken from Refs. [T81IT9] . The addition 
of V^^^ to V^'^ gives a further shift of +0.100 eV. The Is position arising 
from the purely electromagnetic contributions considered so far is —3238.268 
eV. Nearly all the difference between this value and the widely used "total 
electromagnetic binding energy" in Table 1 of Ref. [16] comes from the change 
in the value of — |a^mc, due to the use of the most recent value ^15j of the 
mass He of the 7i~ . 
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There are sizeable contributions to Vr(ls) from the interference of V"^ and 
V^^^ with the hadronic interaction V (see Section 4; as well as the other 
potentials are matrices in a coupled-channel calculation). We determined this 
shift, as well as the width, directly in a subthreshold calculation with mass 
differences included. To this end we evaluated the shift and the width due 
to V, for the two cases V^'' + + V^""* + V and V^^ + V. It turned out 
that, when the direct contributions of V'"''' and V^^* (—3.244 eV and +0.100 eV 
respectively) are accounted for, the shift due to the interference of [V'"^ + V'^^^) 
with V is numerically negligible. The error, estimated by varying the range 
of V, is 0.003 eV. The cancellation of the two interference contributions due 
to V^^ and V^^^* was already noted implicitly in Ref. [16]. The result for the 
width will be given later. 

The shift of the purely electromagnetic Is position given above, due to the 
interference of V with V^'^ only, we call AW. Thus, from the result in the 
preceding paragraph, we have 

W{ls) = -3238.268(8) eV + AW , 

and, using Eq. ([1]), 

AW = -W{3p^ls)+2878.812{8)eV. (2) 

Note that, for practical purposes, the shift AW defined in this way is equal 
to the quantity ei^ defined in Refs. [2|H[T6] . This is because, though we have 
taken account of the short-range effects fully, by including the interference of 
V^'^ with {V""^ + y^^*) in the calculation of the shift and width (rather than as 
a correction to the scattering parameters), the effect is negligible in the case 
of the shift. In the last two equations we have given the uncertainty explicitly. 
It is dominated by the uncertainty in /Xc- 

The hadronic interaction couples the 7r~p channel to the 7r°n and 'yn channels. 
As well as creating the shift AW, it also makes the Is level unstable, with 
a width r. Compared to F, the hadronic and electromagnetic widths of the 
3p level can be neglected numerically. The process then is a three-body decay 
of the well defined 3p level, going via the unstable Is level: 3p {7r^n)'y or 
{'yn)'y. The parentheses indicate that [n^n) and {'jn) are the decay products 
of the intermediate Is level. 

The treatment of such a three-body decay mechanism is given in detail in 
Section 4 of Ref. [20] • Let s be the square of the invariant mass of {ir^n) (or 
(7n)). It turns out that the decay amplitude is proportional to 

(s - + zMoT)-^ . 
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Here Mq = /ic + nip + W{ls), nip being the proton mass. We have neglected 
a possible dependence of F on s. Using the relation between s and the energy 
E of the photon in the rest frame of the decaying 3p state, and noting that 
W{?>p Is) « Mo, one sees that the spectrum of E becomes proportional 
to 



Thus the measurement of this spectrum gives W{'ip Is) and F. 

We now derive the conversion equations which give the values of the elements 
acc and Qco of a, the s-wave scattering matrix for the (7r~p,7r°n,7?7,) system at 
the 7r~p threshold, in terms of W{3p — > Is) and F. We use c, and 7 to denote 
the three coupled channels and qc, qo and k to denote the channel momenta 
in the cm. frame. The matrix a is the value of K at the 7r~p threshold qc = 0. 
K is defined in terms of the unitary matrix S by 



Coiv) 
1 
1 



-Z(S-13)(S + 13) 



Ql/2t-lQl/2 



1 
1 



^ -f3h{7]) 0^ 













(3) 



Equation ([3]) takes account of the low energy behaviour of t induced by 
the long range part of the Coulomb interaction. The quantities appearing 
in Eq. are 



'^gc 0^ 



Q 



go 
A; 



2711] 



exp(27rr7) — 1 



, h{r]) = -\n\r]\ + ^ij{l + ip),{4) 



with P = 2anic , i] = -am^/qc , rric = ficmp/ (fic+nip) . 

The matrix K is real-valued for s > [fie + nip)"^ and can be analytically con- 
tinued into the s-plane. 

The formalism of Ref. [B] will now be used to obtain the conversion formulae. 
Strictly speaking, the formalism relates AW and F to K^c and Kco, evaluated 
at the position of the Is level. For reasons discussed in Section 3, phase shift 
analysis of low energy 7r~p elastic scattering does not give any reliable knowl- 
edge of the range parameters for K^c and Kcq. However, it is sufficient to note 



6 



that the differences between the values of Kcc and Kco at the 7r~p threshold 
and at the position of the Is level are expected to be of order 1 ■ 10~^ fm, 
well below the uncertainties due to the experimental error. We therefore work 
from now on with the matrix a defined above. Following Ref. [6] we introduce 
the matrix A by 



A 



-1 



a 



-1 



^0 0^ 



A; 



(5) 



where go and k are from now on the cm. channel momenta at the tt p 
threshold Qc = 0. Numerically we will use go = 28.0408(96) MeV/c and 
k = 129.40705(31) MeV/c. 

From Eq. we have 

ttcc = ^Acc + go«co'^oo + 2qokacoac^ao^ + Pal^a^^ , (6) 

1 lis 1 

ttco = -|^co|(l + qokal^ + -goOoo - P'^^qi^k^ao^a^^ - -P'^q^kal^) . (7) 



ill 

Equation ([6]) is straightforward; Equation ([7]) uses the result = — ctcoQ'o k~2P~2^ 
where P is the Panofsky ratio, whose value [2T] is 1.546(9). Equations (M7\} 
give only the lowest terms in an expansion. To assess the magnitude of these 
terms we need estimates for Oco, cicy, doo, cto-y, a-yy- The values of these quan- 
tities will be discussed in Section 3. Numerically, we find that 

acc = ^A,,{1 + 4.9 ■ 10-^) , aco = -\A,o\{l + 1.1 • 10"^) . (8) 



The correction terms in Eq. ([8]) are negligible. 



The expression for A 14^ is given in Eq. (12) of Ref. [H] with n = 1, / = 0. 
Keeping only the numerically significant terms we have 

AW 

3?6o(l +pi5Reo[l - i^y] +P2(3?eo)') , (9) 



where 

W{lsr = -^a'm, , ^eo=4ajB , Qe^ = 4qo{l + P-')Ko\V B , (10) 
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B = l/arric being the lowest Bohr radius of the vr p atom. The numerical 
values are = 121.49719(26) MeV, W{1s)p'' = -3234.945(7) eV and B = 
222.56399(48) fm. The quantity pi = 1/2(1 + 7) = 0.78860783, 7 being Euler's 
constant. The expression for p2 is given in Ref. p]. The term involving a 
poorly known effective range parameter is small, and the result p2 = —0.22 is 
sufficiently accurate for our purpose. 

From the experiment [3], Ref. [1] gives W{3p -> Is) = 2885.928(10) eV. Using 
Eq. ([2]), this results in AW = —7.116(13) eV. With this number as the starting 
point for the iteration of Eq. we obtain 

a,c(fm) = 0.0171704(1) (-AW^(eV)) . (11) 

The error in the constant in Eq. (fTTj) comes from the uncertainty in rric (due al- 
most entirely to the error in /ic) and the uncertainty in AW (via the iteration). 
However, it is so small that it is of no consequence and will be disregarded in 
what follows. Using Eq. ([2]), the first conversion formula is 

acc(fm) = (0.0171704) [iy(3p ^ ls)(eV) - 2878.812(8)] . (12) 
Numerically we obtain 

= 0.12218(22) fm . (13) 

We would like to stress that the present experimental results for W{3p — > Is) 
and for F [HIS] are consistent with, but more precise than, the earlier ones 
given in Ref. [2]. The result in Eq. (fT3!) is by far the most accurate number in 
low energy nN physics, and further refinement will not be useful, in view of 
the uncertainty in /ic and even more of the uncertainty in the EM corrections 
to be discussed later. 

We come now to the measured total width F for the decay of the atomic Is 
state into both the vr^n and 'jn channels. The numerical result for the Panofsky 
ratio P quoted earlier is used first to obtain the partial width Fq for decay 
into 7r°n: 

F = (l + p-i)Fo . (14) 

Vacuum polarisation and the finite extension of the charges each have a small 
but significant effect on Fq. The subthreshold calculation discussed earlier in 
connection with the shift shows that V^^^ changes Fq by the factor 0.9948, 
while VP changes it by the factor 1.0049. The total effect of {V^p + 1/"^*) is 
a change by the factor 0.9996. 
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The treatment of decay widths on pages 410 and 411 of Ref. [6] gives 

To = Aa^mlqoal^{l + 2pi5Reo + p?|eo|' + 2p2{'^e^f) , (15) 

keeping only terms which are significant numerically and replacing |AcoP by 
a^o (see Eq. ([8])). Combining Eqs. (1 1411 151) and using numerical values already 
given leads to the second conversion formula 

aeo(fm) = -0.19042(22)(r(eV))^/2 . (16) 



Almost all the error in Eq. (fT6l) arises from the uncertainty in P. An updated 
result for F from the experiment [3] is 0.823(19) eV [5]. This leads via Eq. (fT6!) 

to 

Oeo = -0.1727(20) fm . (17) 



All the error in Eq. (fT7|) comes from the experimental uncertainty on F; the 
error due to the uncertainty of the conversion constant is negligible. The exper- 
iment [3] finally hopes to achieve an accuracy of 1% (an absolute uncertainty 
of 0.008 eV), which will reduce the error on Oco to about 0.0009 fm. Even in 
this case the error due to the conversion constant is only 2 • 10~^ fm; therefore, 
a more accurate measurement of the Panofsky ratio would not improve the 
accuracy of the final value of Ocq. 



3 Calculation of electromagnetic corrections 



To calculate the electromagnetic corrections to a^c and Oco due to the point- 
charge Coulomb interaction, the external mass differences and the presence 
of the 7n channel, we set up a potential model for the physical three-channel 
situation at the 7r~p threshold. The relativised Schrodinger equation (RSE) 
for the s-wave has the form 

[S/drH^ + Q2 - 2i5i(V + VP^)]u(r) = , (18) 



where, since we are working at the vr p threshold gc = 0, 



Q 



'^O 0^ 

go 
yO /cy 
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with go and k having the same meaning as in Eq. ([5]). The diagonal matrix m 
of modified reduced masses is 



^m, 0^ 



m 



mo 
m, 



(19) 



For a two-body channel with masses mi and m2, the modified reduced mass 
is 



m 



W"^ — ml — ml 
2W 



Thus, for W = He + rn.^ 



/ic + mp 



mc, mo 



(/ic + mpY - fil-m: 



2(/ic + mp 



,m^ 



(/ic + mp)^ - ml 
2(/ic + mp) 



(20) 



The reasons for the inclusion of these relativistic modifications of the reduced 
masses for the 7r~p and ir^n channels were discussed in detail in Section 3 of 
Ref . [22] , for the two-channel potential model calculation of the EM corrections 
for the analysis of low energy n~p scattering data. The modification is essential 
for the 7n channel, where the standard reduced mass is zero, whereas m^ in 
Eq. (!20l) is in fact just k. 



The 3 X 3 matrix V^^ in Eq. ([T8]) is 



V 



pc 








where V^'^ is the attractive point-charge Coulomb potential. We recall that 
the effects of vacuum polarisation and of the finite extension of the charges of 
TT~ and p have already been accounted for. We write the matrix V in the form 



V 



^ Vcc VcO Vcj ^ 
VcO Voo Vqj 

y Vc-y Vq^ V^^ j 



and define the submatrix 'V ^ - 




(21) 
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We shall call the EM modified hadronic potential because, as explained 
in Section 1, it contains residual EM effects. There is clear evidence, to be 
discussed in Section 5, that violates isospin invariance, though only by a 
few percent. However, we have only the numbers a^c and Qco available exper- 
imentally (ooo is inaccessible) and therefore need to express Yh in terms of 
two parameters. Hence we made the approximation of treating V/i as isospin 
invariant: 



We made the same approximation in Ref. [22] in order to calculate EM correc- 
tions for the analysis of low energy 7r~p scattering data. Since this potential 
is used solely for the calculation of the corrections, this approximation should 
have very little effect on the results. For V1/2 and V3/2 we assumed the radial 
dependence of the s-wave hadronic potentials used in Refs. p^|22] and left 
their strengths as parameters to be adjusted. The leading term in each po- 
tential is of gaussian form, with a range parameter which is fixed at 1 fm (a 
realistic value for the hadronic interaction), and two small correction terms 
are added. 

The three potentials 14^, Vo-y and V^^ in Eq. (!2Tll also need to have some as- 
sumed radial dependence, with their strengths left as parameters to be ad- 
justed. A study of the analytic structure of the relevant amplitudes shows 
that the range parameter 1 fm is also appropriate for these three potentials. 
For convenience we took the radial dependence of to be the same as that 
of Ko, and the radial dependence of V07 and V^^ to be the same as that of 
Voo- In the calculation there were then five potential strengths to be adjusted, 
with five experimentally determined numbers used to fix them. 

In addition to the values of Occ and given in Eqs. fll3|17p . we used experi- 
mentally determined values of Oc^, ao^ and a^^. From the values of and P, 
we obtain = 0.0647 fm. The value of oq^ can be obtained only from an 
analysis of experimental data near threshold for the photoproduction process 
7n — i> vr^n. The results given in Ref. [23] lead to ao-y = — 0.0028 fm, but no 
error is given; we think that 0.0006 fm is a realistic uncertainty. The value of 
a-y^ comes from the analysis of data on Compton scattering off neutrons. From 
results in Ref. |24j we have extracted a^^ = 0.0011 fm, with an uncertainty 
of around 20%. The five potential strengths were then varied until they re- 
produced the five numbers just discussed, thus fixing V for the three-channel 
calculation. 

We now give some details of the calculation of a from the RSE (ITSl) . For 
the two-channel case (omitting the 'yn channel) the calculation is given fully 
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in Section 3 of Ref. [25]; it can easily be generalised to the three- channel 
situation. The first step is to form the matrix P(r) by putting side by side, 
as its three columns, three linearly independent regular solutions u^{r),j = 
1, 2, 3, of Eq. (ITHI) . Let i? be a distance beyond which V is negligible. We then 
define the real symmetric matrix d by 

d = Rm-^/^P'{R)[P{R)]-^m^/^ , 



where the prime denotes differentiation with respect to r and the matrix ele- 
ments of m are given in Eq. (120|) . Note that our present notation is different 
from that of Ref. [25] . The matrix K in Eq. (64) of that reference is now t, 
defined just before Eq. (|3]), and A in Eq. (65) is now K as defined in Eq. ([3]). 
Equation (66) of Ref. ^], generalised to the three-channel case and taken at 
threshold, gives 



R 



-1. 



X 



- 4oPo ^ sin po cos po - c^oo Po ^ sin po - (io^Po ^ sin po 

—dc-^fP~^ sin p^ —do^p'^ sin p^ cos p^ — d^.yp~^ sin p^ 
/ 

-'^o + dccXo dcoXo 

dco cos Po Po sin po + doo cos po d, 



dc-yXo 



07 cos Po 



v 



dcj cos pj 



d 



COS Py 



Py sin Py + dyy COS Py ^ 



:22) 



where po = qoR and py = kR. The auxiliary Coulomb functions $o, '^'o, -^o 
and Xq (discussed at length in Section 2 of Ref. ^25]) have at threshold the 
arguments {—PR; 0), which we have suppressed. In terms of Bessel functions 
they are, with x = 2^/J3R, 

$0 = 2x~'Ji{x), = Mx), Xo = -^TTxY^ix), X* = ^x^Yoix). (23) 



Having given the details of the modelling of the three-channel situation, we 
now proceed to calculate the corrections which remove from Occ and Oco the 
effects of the point charge Coulomb potential, of the mass differences pc — Po 
and rrin — mp, and of the presence of the •yn channel. To obtain them, we need 
to model the two-channel situation with these three effects removed, using the 
RSE at threshold 



[d^/dr^l2 - 2m^Vh]u{r) = 



(24) 



The appearance of in Eq. (12^ implies the assignment of the mass Pc to 
both 7T~ and 7r°, and the mass nip to both p and n. This choice is made in 
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most studies of the low energy nN system and was discussed fully in Section 
5.2 of Ref. [9]. Thus the scattering parameters (ice and (ico which result from 
Eq. (12^ still contain residual EM contributions. 

By integrating Eq. fl24p we obtain the 2x2 matrix a; the numerical values for 
the relevant matrix elements are 

= 0.12140 fm , = -0.1757 fm . (25) 
We now define the stage 1 EM corrections 

A^flcc = dec - Cicc , A^OcO = ttcO - CLcO ■ (26) 

From Eqs. (DSD, (HI]), (I2SD and ([2S]), we have 

A^Occ = +0.00078 fm , A^a^o = +0.0030 fm . (27) 

These are the essential results from the present calculation. 

It is instructive to look at the contributions of the three effects to the correc- 
tions A^acc and A^Ocq. For A^acc we have 

Aaccipc) = -0.00214 fm,Aacc(md) = -0.00071 fm , Aacc(7n) = +0.00363 fm. (28) 

The smallness of the total correction is thus seen to be the result of a near 
cancellation between the effect of the jn channel, which is large, and the 
combined effect of V^'^ and the mass differences. One sees how important it is 
to estimate the effect of the jn channel as reliably as possible. For A^a^o the 
situation is quite different. We now have 

Aaco(pc) = -0.0002 fm,Aaeo(md) = +0.0029 fm , Aaco(7n) = +0.0003 fm. (29) 

The effect of the 7n channel is now very small, and the correction is due almost 
entirely to the mass differences. 

Uncertainties in the corrections given in Eq. fl27j) arise from: (a) the shapes 
assumed for the components of V/^; (b) the errors of the five experimentally de- 
termined quantities which are fitted; (c) the approximation of using an isospin 
invariant form for V^. The quantity Aaccipc) has a significant uncertainty due 
to lack of knowledge of the shapes of the components of V/j. This was esti- 
mated as in Refs. p^1[22] by varying the range parameter between 0.8 fm and 
1.2 fm and simultaneously adjusting the strengths of the potentials so that 
the five scattering parameters discussed earlier in this section remain constant. 
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The quantity Aacd'^n) suffers from the same uncertainty, and has in addition 
a significant error due to the uncertainty and a^^. The uncertainty 

in Aaco(^«(i), the dominant part of A^Oco, comes from uncertainties in V/j (the 
shapes of V1/2 and V3/2 and the assumption of isospin invariance). 

Unfortunately, in the situation under consideration one cannot achieve the 
precision associated with the errors quoted in the experimental results. The 
variations described in the previous paragraph can lead to rather rough esti- 
mates only. It is therefore wise to err on the side of caution and increase these 
estimates accordingly. We think that, on the basis of the estimates we have 
made, the errors on each of the components of the corrections could be as high 
as 20%. This leads to the following final values for the corrections 

A^Occ = +0.0008(8) fm , A^aco = +0.0030(6) fm . (30) 

Combining Eqs. f|T3|) . f|T7j) and fl30l) . we obtain the final results 

Occ = 0.1214(8) fm , 5co = -0.1757(21) fm . (31) 

(We have given results for the s-wave scattering parameters in fm because it 
is a precisely defined unit. In units which are in common use, the final results 
in Eq. ^ are Sec = 0.0859(6)/^-^ and a^o = -0.1243(15)/!"^) 

It is important to note that the error in ficc arises almost entirely from the 
uncertainty in A^a^c, so that an improvement in the accuracy of the measure- 
ment of W{?)p Is) will not reduce the error in a^. On the other hand, 
almost all the error in 5^0 comes from the experimental uncertainty in F. If 
the error in the measurement of F can be reduced to 1%, the uncertainty in 
Sco in Eq. (13T|) will be reduced from (21) to (11). 

In the discussion following Eq. (35) of Ref. [9j|, we emphasised the very large 
difference between the values of (ice which come from pionic H and from our 
PSA of the combined vr+p and 7r~p elastic scattering databases for < 
100 MeV. Converting from the unit to fm, the latter gives the value 
Sec = 0.1127(16) fm. The difference between this value and the value 0.1214(8) 
fm given in Eq. (131 p is thus 0.0087(18) fm. This large difference is a subtle 
consequence of the violation of isospin invariance in the electromagnetically 
modified hadronic interaction, which is convincingly demonstrated in Section 7 
of Ref. |9] and in Section 5 of this paper. The phase shift 5i/2 used in Ref. p] to 
analyse 7r~p elastic scattering data is an artificial construct, which enables the 
analysis to be done in a framework of formal isospin invariance. When isospin 
invariance is violated, 81/2 goes to a nonzero value at threshold, so that the 
scattering length ai/2 given in Ref. ^ is an artefact of the parameterisation of 
(5i/2, which forces it to behave like qc near threshold, and the extrapolation to 



14 



threshold (using this artificial parameterisation) of values of (5i/2 obtained from 
data at momenta qc ~ 80 MeV/c upwards. Therefore the value of acc extracted 
from pionic H cannot be compared directly with the value of (2ai/2 + 5+)/3 
obtained in the PSA of Ref. [H]. 



4 Comparison with other calculations of the electromagnetic cor- 
rections 

As mentioned already in the introduction, in our method of calculating the 
EM corrections to the pionic H data we use hadronic potentials which are of 
the same form as those of Refs. [T01I22] . This implies that the EM corrections 
which were applied to the low energy vr^p scattering data in our PSA of Ref. |9] 
are compatible with the ones obtained in the present paper for the pionic H 
case. For the PSA, it was sufficient, in analysing data from ~ 80 MeV/c 
upwards, to take the 7n channel into account by using the perturbative result 
of Ref. [26]. It has practically no effect on the results of the PSA. However, for 
our calculation at threshold, the effect of the 7n channel turns out to be very 
significant in the case of the correction to Occ, thus justifying the full three- 
channel calculation we have made. In these respects, our calculations of the 
EM corrections to the pionic H data differ considerably from the older work 
of Refs. [TBPT] . both of which use primitive forms of the hadronic potentials 
which reproduce only the scattering lengths, in contrast to our forms which 
also reproduce the phase shifts for T^^ < 100 MeV. Furthermore, Ref. [27] is 
a nonrelativistic calculation which uses much older values of the scattering 
lengths. In Ref. p2j the 7n channel is taken into account in an unsatisfactory 
way, by using the Panofsky ratio to alter the input value of aco in what remains 
a two-channel calculation, whereas we have included this channel in the most 
realistic way possible in a three-channel calculation. In addition, Refs. p^][27] 
do not take account of the accurate relations given in Ref. [6], between AW 
and ttcc and between Pq and Oco- This means that our corrections and those 
of Ref. [in] differ by the factors in the parentheses of Eqs. (piT^ . For these 
reasons, we consider that our calculations supersede the earlier potential model 
calculations in Refs. [T6|l27j . 

A quite different method of calculating the EM corrections is used in Ref. |13j . 
In a nonrelativistic approach using Coulomb wave functions, with a short- 
range hadronic interaction and extended charge distributions of vr" and p, the 
authors identify four components of the corrections. Two are the contributions 
from the interference of V^''^ and V^^^ with the hadronic potential V; the other 
two are a renormalisation term and a gauge term. The interference terms 
are obtained using an argument based on the change at the origin in the 
nonrelativistic Is wave function arising from V^'^, when V^^ and V"^^* are 
added to V^'^. This is a purely EM calculation, and the results are therefore 



15 



claimed to be independent of V. We can see no theoretical justification for 
this procedure of artificially introducing the wave function at the origin as 
a replacement for the lowest Bohr energy in the DGBT formulae, and then 
calculating corrections via changes in this wave function. Moreover, while the 
shift in the Is level position due to the interference of V^^ with V happens to 
be correct (—0.0034 eV), that for V^^^* comes out to be +0.0060 eV, whereas 
the correct result is +0.0034 eV, which gives the cancellation reported in 
Section 2 and in Ref. [16]. The subthreshold calculation discussed in Section 
2 shows that these contributions separately have a significant variation with 
the range of the hadronic potential, and we have given the error which results 
from this for the sum of the contributions. Ref. [13] also incorrectly says that 
the results given there for the interference contributions agree with those given 
in Ref. 116]. 

While Ref. [13] has a discussion of the multichannel situation in Subsection 
2.2, it does not explicitly give the terms which appear in parentheses in our 
Eqs. ( MTSj) . nor are these (admittedly small) terms taken into account in 
their corrections. Further, we cannot see any corrections in Ref. [13] which 
correspond to the corrections due to external mass differences appearing in 

Eqs. mm. 

That leaves the renormalisation and gauge terms in Ref. [TB]. Each of these 
terms involves an integral over the extended charge distributions, so it appears 
that some of the effect of extended charges is contained in these terms. In our 
approach, however, V^^^* is fully accounted for by its direct contribution to 
the position of the Is level and the interference of V^^^ with V. It also does 
not seem possible to see where the interference of V^'^ with V, which gives 
the corrections designated (pc) in Eqs. ( |28|29ll . appears in the corrections of 
Ref. . So here we deal with two very different approaches to the calculation 
of the EM corrections, with no discernible way of comparing them. 

Further, the gauge terms in Ref. \13\, for the corrections to and Oco, are 
much more uncertain than is claimed there. They are proportional to the range 
parameters b'^^ and b^Q respectively, which appear in the low energy expansions 
of Kcc and Kcq. Ref. [13] uses for these quantities values given in Ref. |[28| . 
taken from a PSA of ttN scattering data obtained before 1983. Practically 
none of these data were used in our recent PSA [9]. The numbers in Ref. [28] 
are no longer of any value. Moreover, as we explained at the end of Section 
3, because of the violation of isospin invariance in the electromagnetically 
modified hadronic interaction, the PSA of Ref. [H] does not lead to a reliable 
estimate of dec- For b^c the situation is even more uncertain. The value obtained 
from the PSA cannot be accepted as a reliable determination; it is merely a 
rough guess based on phase shifts for > 80 MeV/c. The same holds for b^Q. 
It is not possible at present to obtain reliable values of the range parameters 
from a PSA. As a result, the values of the EM corrections proposed in Ref. [13] 
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are quite uncertain. 



In our calculations, the range parameters are replaced by potentials which, as 
explained in Refs. [TU|I^ . are constructed to reproduce the phase shifts for 
Qc ^ 80 MeV/c. The lack of experimental data at very low energies increases 
the uncertainty in the potentials; it allows a wider interval of values of the 
gaussian range parameter a which yield acceptable fits than would otherwise 
be the case. We have taken this uncertainty in a into account in the estimates 
of the errors in our corrections given in Eq. fISUl) . 

In Ref. [29] the dispersive EM correction to Occ which arises from the radiative 
processes 7r~p — > '~^X{X = n, A) is calculated in the heavy baryon limit. Table 
1 of Ref. [22] gives the correction Aacd'yn) = 0.00424(13), while our result 
(Eq. (12511 ) is 0.00363(72) fm. These two results, which have been obtained in 
completely different ways, agree quite well. The case X = A is very different. 
In the potential model, the effect of a 7A intermediate state is a residual EM 
contribution to V, and has to be calculated in another way. The result given 
in Table 1 of Ref. [29j is Aacd'jA) = +0.00332 fm. There is no other estimate 
of this correction. One should note that diagrams with a 'ynN intermediate 
state in principle appear in the calculation of Ref. [30] , which we shall discuss 
in a moment, but at a higher order in the perturbation expansion. An effect 
due to 7A as large as that given in Ref. [29] would imply that higher orders 
of the expansion give contributions which are comparable in magnitude to the 
one-loop result given in Ref. [SO] . 

Gasser et al. [30] have made a calculation of the EM correction to Occ us- 
ing baryon chiral perturbation theory (ChPT) in the infrared regularisation. 
The calculation was carried out up to diagrams containing one loop; it is a 
calculation of the complete EM correction to to this order. An unam- 
biguous splitting of this correction into stage 1 and stage 2 components is a 
very difficult task, and may not even be possible (for example, due to the 
ultraviolet divergences). Nevertheless we can say that some parts of the cor- 
rection correspond to the effects calculated (to all orders) with our potential 
model. In addition, there are much more important parts of the correction 
of Ref. [30|, which are not taken into account in our calculation. They give a 
large negative correction which swamps the small correction A^acc in Eq. fl30|) . 
Unfortunately, this correction is very uncertain, due to a lack of knowledge of 
the low energy constant /i. 

For the width, there is a calculation using baryon ChPT reported in Ref. |31j . 
but only to lowest order (no loops). It contains contributions which are not 
accounted for in the potential model calculation; the total contribution is 
positive, but very small compared with the correction A^Oco in Eq. (130]) . 
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5 Violation of isospin invar iance 

If the electromagnetically modified hadronic interaction were to satisfy isospin 
invariance, we would have the equality 

ttcc + V^CLco = a+ , 



where a+ is the s-wave scattering length for elastic vr+p scattering. The value 
of a+ needs to be determined by means of a PSA of low energy iT~^p scattering 
data, using EM corrections calculated in a manner consistent with the cal- 
culation of the corrections in Eq. (!30i) . Our analysis [9] of the n'^p scattering 
data for < 100 MeV gives 

S+ = -0.1062(55) fm . (32) 



When the analysis is made of a database which combines ir'^p and 7r~p elastic 
scattering data up to 100 MeV, using a more restrictive model for the fitting, 
the value of 5+ is practically unaltered, but the error decreases from (55) to 
(23). For the purpose of testing for isospin invariance, the larger uncertainty 
in Eq. fl32|) should be used. 

From the results in Eq. (|3T|) we have 

ttcc + V2aco = -0.1271(31) fm . (33) 



Almost all the error in Eq. (155]) arises from the uncertainty in a^o- Reduc- 
ing the error in F to 1% would reduce the uncertainty in Eq. fl33|) to (17). 
The difference between the numbers in Eqs. f l32|33l) is 0.0209(63) fm, which 
amounts to 3.3 standard deviations. This is evidence for a significant violation 
of isospin invarance at the electromagnetically modified hadronic level. If the 
experiment p] achieves its hoped for accuracy for the measurement of F, the 
error in the difference will change from (63) to (58). If the result for F remains 
approximately the same, the difference will become even more significant. 

This result strengthens the evidence for the violation of isospin invariance 
at the electromagnetically modified hadronic level given in Ref. [9], which 
comes from the inability of the phase shifts obtained from the analysis of the 
combined n'^p elastic scattering data to reproduce the data on the charge 
exchange reaction 7i~p — > n^n. 
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